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Abstract 

The phase structure of a non-isotropic non-Abelian SU(3) lattice gauge model 
at finite temperature is investigated to the third order in the variational-cumulant 
expansion (VCE) approach. The layer phase exists in this model in the cases of 
dimensions D = 4, D = 5 (d = D — 1). 
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1 Introduction 

X 

Up to now , though all experiments do not deny that the universe people live in a (3+1)- 
dimensional space-time, we need not to believe absolutely that our universe must be so|IJ. 
Perhaps our universe is only a domain of a higher dimensional space-time, which we do not 
realize at this moment. This thought is inspired by the quickly development of string the- 
ory recently. Consequently, there must be critical differences in physics between our (3+1)- 
dimensional universe inside and outside. New phenomena need to be investigated to some 
extent. 

An interesting non-isotropic Abelian lattice gauge model on the existence of a high di- 
mensional space-time has been proposed in Ref. || by one of us, Fu, and Nielsen, where the 
coupling constants of gauge interactions depend on the spatial-time directions. A new phase, 
(i-dimensional layer phase, has been found in a D-dimensional {D > d) non-isotropic U(l) 
lattice gauge model at zero temperature by using mean field method and strong coupling cor- 
rections. The results have been confirmed by Monte-Carlo (MC) simulation (3| and by other 
method, e.g., by using the variational-cumulant expansion approachQ. Further investigation 
show that the layer phase exists in the non-isotropic Z(2)[^], non-Abelian SU(2)f|] (confirmed 
by MC simulation HI) and SU(3)|| lattice gauge model, also exists in the non-isotropic U(l) 
lattice gauge model with bosons and fermions[^, ||] at zero temperature. Furthermore, recent 



works confirm the exists of layer phase in non-isotropic U(l) lattice gauge model in Ref[10] and 
show that the layer physe exists in the non-isotropic U(l) lattice gauge model with Higgs [|l~l]] . 
The phase structure of non-isotropic U(1)|I^| and SU(2)fl3]l lattice gauge model at finite 



temperature has been studied with the variational-cumulant expansion model, the layer phase 



is also found which is different from the zero temperature case. In this paper, we shall use 
the variational-cumulant expansion approach [ 14 1 to investigate the phase structure of the non- 
isotropic non-Abelian SU(3) lattice gauge model at finite temperature. 



2 The non-isotropic SU(3) theory using the VCE at 
finite temperature 

If we suppose that there is a D-dimensional space-time described by SU(3) gauge symmetry 
with the coupling constants depending on the spatial-time directions. Then the Wilson action 
of this non-isotropic SU(3) lattice pure gauge model is: 

5 = E§7 E ^p + C)=E| E Tr(U p ^ + U^), (1) 

where U Palxv is an ordered product of Ui €SU(3) around a plaquette, \i and v are the space-time 
directions, N = 3 is the rank of SU(3). Since the key point of the finite temperature lattice 
gauge field theory is to establish an appropriate partition function with a periodic boundary 
condition of the field in the time- like direction(s), the time- like link(s) and space-like links must 
be treated differently, in the non-isotropic SU(3) lattice gauge model at finite temperature, we 
suppose that there is only a one-dimension time-like direction which is defined along the dth. 
direction in the D-dimensional non-isotropic space-time lattice. Furthermore, we suppose 
there is a d-dimensional subspace which are d dimensional space-time including the time-like 
direction in the Z?-dimensional space-time lattice. In the d dimensional subspace the coupling 
coefficients are assumed identical and called /3, the remaining [D — d) coupling coefficients are 
also taken to be identical and we call them 7. The reason we choose two different coupling 
constants (3 and 7 are that in one side our familiar universe is an isotropic system, and on the 
other side we mainly attempt to show the differences between the non-isotropic gauge system 
and our familiar isotropic gauge system. Having these supposition, the Wilson action of the 
non-isotropic SU(3) lattice gauge model at finite temperature can be written as: 



2N 



£ T T (U Pa + U+)+ £ Tr(U PT + U+] 



+ 
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2N 



E Tr(U p , a + U+)+ £ Tr(U p ,+U+) 



(2) 



here one may consider the first term as an "inside layer" and the second term represents an 
interaction between different layers. Where fx, < u = 1, 2, • • • , d — 1), //(// = 1, 2, • • • , d — 
1, d + 1, • • • , D) and v'{y' = d + 1, • • • , D) are space-like directions, p a and p a i are the space-like 
plaquettes, and p T , p T > are the time-like plaquettes. The summations {• • •} are taken over all 
such plaquettes. 

The corresponding partition function of the system is: 



[dU]e i 



(3) 



In general, it is very difficult to calculate the partition function of a certain system. According 
to the variational-cumulant expansion method, if one introduce a trial action, the calculation 
of partition function will become simple. We introduce a trial action: 

S = £ Tr(J+U a + J a U+) + £ Tr(J+U a , + J a ,U+) 

a a i 

+ Y,Tr(J+U T + J T U+), (4) 

T 

J CT , J T and J a > are real variational parameters corresponding to the space-like link variable 
U£, the time-like variable U^~ in the cf-dimensional subspace, and the spacedike variable U^i 
in the remaining (D-d) directions, respectively. For the convenience of calculation, we let 
J a = Jt = J and (Ja)ij = Jij = J$ij and (J a ')ij = Ja'^ij- The corresponding trial partition 
function of the virtual system can be calculated by VCE^j: 

Z(J, J' a ) = J [dU]e s ° (5) 

Then the partition function of a genuine system can be achieved via the trial system as: 

Z = J [dU]e s = J [dU]e So e s - So = Z(J a , J ff ,) (e S ~ S °) o , (6) 

where 

expresses statistics average in the trial system. 

Now, let us introduce the order parameter, Polyakov line, which can be used as a criterion 
of phase transition of a given system: 

( L ) = (^ Tr (E 17 ^))) = \ J ^ LeS 

/ OC | CO -i 

[dU]Le s - s °e s « = -i(L(S - S ) n )c = £ 3 <A»)o (7) 

n=0 n=0 

where N T = 2 is the periodic condition, (• • -) c indicates the cumulant average in the trial system 
which can be expanded in the statistical averages and the lower order cumulant averages 1 14]. 
For example 

(Lq)c = {L)o, 

(Li)c = (L(S - 5 )) c = {L(S - 5 ))o - (£}o<5 - S )o, 
(L 2 ) c = (L(S - S ) 2 )c = (L(S - S ) 2 )o ~ 2(5 - S )o(L(S - S )) c - (L) ((S - S ) 2 )o, 

(8) 

In this work, we calculate the order parameter Polyakov line to 3rd order (n=2). 



3 Variational treatment 



In principle, we can calculate the Polyakov line to any desired orders using the VCE method 
and obtain the exact Polyakov line of the real system without any dependence on the variational 
parameters J a and J a i. In practice, we can only expand (L) to some limited terms and the 
expansion Eq.(8) must be truncated at some order n(in our work, the (L) is expanded to 
n = 2). Thus, the (L) turns out to depend on J a and J a i. Then the choice of the parameter J a 
and J a i will affect the rate of convergence of the expansion. We use the convexity inequality 
to obtain the best values of J a and J CT /. First, we introduce the equation of free energy, 

Z = Z(J a , J a ,) (e s - s °) Q ^ e~ w = e~ w ° (e s ~ s °) o 

where the W and Wo = —ln[Z(J a , J a >)o] represent the free energy of the real system and trial 
system. Then we get the standard convexity inequality equation is 

W = -ln[Z] = -ln[Z(J a , J a ,) (e 5 " 5 °) o ] < -ln[Z(J a , J a ,) } - (S - S )o = W eff , 

where W e ff is the upper bound of the free energy. 

So, the variational parameters J a and J' a are determined by minimizing the W e ff. 

which lead to 

where Z (a) and Z (a r ) are the single link group integral, for example, the definition of Z (a) 
is 

Z (a) ee Z (J a ) = J [dU]e T ^ u++UJ t\ 
Z' Q (a) and Z (a') are the corresponding differential of Zq{u) and Z${<j'). 
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4 Result 

Eq.(lO) have solutions in the following forms: 

(1) J a = (Z {J a ) = 0, Z^J a ) = 0), J a , = (Z (J a ,) = 0, Z {J a ,) = 0), 

(2) J a > (Z (J a ) > 0, Z' {J a ) > 0), J a i = (Z (J a ,) = 0, Z {J a ,) = 0), 

(3) J a > (Z (J a ) > 0, Z (J a ) > 0), J u' > (Z (J a ,) > 0, Z (J a ,) > 0). 
Corresponding to the three kinds of solutions there are three kinds of phases: 

(a) The high-temperature phase corresponding to solution (1). The order parameter is: 

(L) = 0. 

(6) The layer phase corresponding to solution (2). The order parameter is: 



(L) = (L((3)}. 



(c) The low temperature solution corresponding to solution (3). The order parameter is: 



(L) = (L%i)). 

Here (L) is calculated to the terms of 1 and 7 2 and the calculation is very tedious, we will 
not show it here. 

Base on the VCE approach to 3rd order, we obtained the layer phase diagrams of the 
non-isotropic SU(3) lattice gauge model at finite temperatures shown in Fig.l and Fig. 2. As 
indicated by the lattice gauge field theory of zero temperature §|, the layer phase of d = 
D — 1(D < 5) is unstable, but stable when d = 5(D = 6) at zero temperature. While in our 
calculation of finite temperature the layer phases exist and are stable when d = D — 1(D = 4, 5) 
for a fixed N T value (here N T = 2), which corresponds to a not very low temperature, as 
iV T — * oo for zero temperature. It implies that there must be a transition temperature at 
which the layer phase turns to be unstable. Unfortunately, for the case of larger iV T value, 
the calculations with VCE are too complicated and difficult to carry out. The triple points 
of phase diagram in Fig.l and Fig. 2 are located at j3 c = 9.2, 7 C = 3.9 and C = 6.2, 7 C = 3.1, 
respectively. 

We do not calculate the cases when D > 6, but by our experience, the layer phase will exist 
in those cases. For conclusion, using VCE approach, we get information about the layer phase: 
the layer phase exist in the non-isotropic non-Abelian SU(3) gauge model, which conform the 
former imagination |4j] that the layer phase not only exist in the non-isotropic U(l) and SU(2) 
gauge model, but also in other non-isotropic gauge models. But in this work we can only 
obtain that the layer phase exists at a very high temperature in the non-isotropic SU(3) gauge 
model, we can not give the exact temperature at which there is phase transition and the layer 
phase disappear. So, further investigation is needed. 
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Figure 1: Phase diagram for the non-isotropic SU(3) gauge model when D = 4 and d = 3. 
The upper-right area shows the low-temperature phase, the upper-left area shows the layer 
phase and the lower area shows the high-temperature phase. 



14 I i i i i | i i i i | i i i i | i i i i | i i i i 

P : 




Figure 2: Phase diagram for the non-isotropic SU(3) gauge model when D = 5 and d = 4. 
The upper-right area shows the low-temperature phase, the upper-left area shows the layer 
phase and the lower area shows the high-temperature phase. 
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